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CHAPTER  1 


INTRODUCTION 


1 ..L  Motivation  for  the  Use  of  Microstrjp  Antennas 

In  recent  years  there  has  been  an  increased  interest  in  microstrip 
antennas,  primarily  due  to  their  surface  conformability , compact 
structure  and  ease  of  construction.  A general  microstrip  antenna  can  be 
pictured  as  in  Figure  1.1.  The  antenna  element  is  separated  from  the 
ground  plane  by  a suitable  low-loss  dielectric  medium.  Since  these 
antennas  have  such  a low  profile,  they  are  suitable  for  high  velocity 
vehicles  where  wind  drag  must  be  kept  to  a minimum. 

Experimentally,  these  antennas  have  been  found  to  have  a relatively 
narrow  bandwidth.  Yet  because  these  antennas  are  so  thin,  there  is 
still  active  interest  in  finding  a simple  analytic  model  to  describe  and 
predict  their  behavior. 

A major  use  for  the  microstrip  antenna  is  as  an  element  in  an  array. 
Experimental  results  indicate  that  no  "magic"  element  shape  exists  that 
will  give  a substantially  wider  impedance  bandwidth  than  any  other 
shape.  For  this  reason,  it  would  seem  more  important  to  concentrate  on 
a short  list  of  simple  geometries,  rather  than  a long  list  of  random 
geometries.  Of  course,  the  short  list  must  include  elements  whose 
radiation  patterns  are  sufficiently  varied,  so  as  to  allow  a broad  range 
of  array  design.  Such  a list  could  include  the  disk,  half  disk, 
annulus,  rectangle,  square  and  triangle.  A simple  and  accurate  theory 
has  been  developed  for  use  with  these  geometries.  In  Chapter  2,  the 
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disk  element  is  analyzed,  and  the  method  of  analysis  can  also  be  applied 
to  other  simple  geometries. 


L.2  Literature  Review 

In  the  literature,  some  scattered  and  limited  attempts  have  been 
made  to  analyze  the  behavior  of  a few  microstrip  elements.  Some  of 
these  attempts  have  incorporated  doubtful  and  even  invalid  assumptions. 

Munson's  analysis  of  the  rectangular  element  is  one  example,  [1]. 
It  was  claimed  that  the  input  impedance  could  be  modelled  as  two  slot 
antennas  separated  by  a low  impedance  transmission  line.  The  obvious 
oversight  is  that  this  uses  superposition  of  the  power  patterns  of  the 
two  slots,  thus  neglecting  any  mutual  coupling.  A similar  approach  was 
used  by  Derneryd  [2].  No  attempt  was  made  to  compare  calculated  and 
measured  input  impedance  results  in  either  [1]  or  [2],  and  their 
analysis  is  very  limited  since  the  feed  point  must  remain  on  that  edge 
of  the  rectangle  where  the  E field  is  constant.  Finally,  the  method  can 
not  be  used  for  any  other  element  shape. 

Another  method  [3]  used  a grid  of  linear  and  V shaped  dipoles  to 
approximate  the  radiating  element.  This  method  shows  an  agreement 
between  calculated  and  measured  input  impedances  that  is  quite  inferior 
to  that  obtained  by  the  simple  analysis  in  Chapter  2.  Even  so,  the 
method  was  inefficient  with  respect  to  the  computer  time  required.  The 
analysis  of  Chapter  2 shows  excellent  results,  and  requires  orders  of 
magnitude  less  computer  time  than  that  of  [3].  As  an  example,  the 
resonant  radiation  pattern  and  the  input  impedance  for  15  frequencies 
near  resonance  were  calculated  in  under  2 seconds  for  a disk  antenna. 


CHAPTER  2 
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ANALYSIS  OF  THE  MICROSTRIP  DISK 


2,1  Th<?  Simple  Cavity  Model 

The  method  used  in  analyzing  the  disk  is  somewhat  different  from 
methods  found  in  the  literature,  and  is  quite  simple.  Analytic  results 
obtained  using  this  method  agree  very  well  with  measured  data,  and  it 
appears  that  validity  has  not  been  sacrificed  for  simplicity. 

Power  is  fed  to  the  antenna  by  either  a coaxial  line  piercing  the 
ground  plane,  or  by  a microstrip  transmission  line  at  the  edge,  as  shown 
in  Figure  2.1.  The  antenna  is  modelled  as  a closed  cavity,  bounded  on 
top  and  bottom  by  perfect  electric  conducting  (PEC)  planes,  and  around 
the  perimeter  by  a perfect  magnetic  conducting  (PMC)  ribbc  r,  as  in  Figure 
2.2.  Henceforth,  the  disk  will  be  oriented  with  the  z axis  as  shown. 

The  PMC  wall  approximation  would  seem  plausible  with  respect  to  the 
fields  between  the  conducting  planes  since  the  radially  directed  surface 
current  (on  the  underside  of  the  top  PEC  plane)  must  approach  zero  at 
the  edge.  Of  course,  this  can  not  be  assumed  for  the  external  fields 
since  such  a boundary  would  not  allow  radiation.  However,  the  radiated 
power  can  be  approximated  once  the  internal  fields  are  known  at  the 
edge.  The  source  of  the  radiation  is  taken  to  be  the  edge  distribution 
of  the  internal  fields.  Calculating  the  radiated  power  in  this  manner 
requires  the  assumption  that  the  internal  field  structure  is  relatively 
independent  of  the  radiated  power.  The  measured  quality  factor  of  these 
microstrip  antennas  is  high  so  that  this  assumption  is  plausible. 
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Since  the  antenna  is  relatively  thin,  the  internal  fields  are  taken 
to  be  TEM.  Indeed,  for  the  frequencies  of  interest,  the  TE  and  TM  modes 
of  all  orders  are  well  below  cutoff  and  the  electric  field  is 
perpendicular  to  the  disk  surface,  with  no  divergence. 

When  excited,  the  internal  field  structure  of  this  simple  cavity 
model  can  be  found  using  either  mode  matching  or  a modal  expansion.  It 
is  then  straightforward  to  calculate  the  stored  energy,  radiated  power, 
internal  copper  loss  and  the  driving  point  voltage.  When  all  of  these 
parameters  have  been  determined,  the  input  impedance  and  radiation 
pattern  can  be  computed.  It  is  because  the  geometries  are  simple  and 
fit  nicely  into  some  coordinate  system  that  mode  matching  and  modal 
expansion  can  be  used  here.  This  is  one  motivating  factor  for 
concentrating  on  a short  list  of  simple  element  shapes. 


2,2  The  Internal  Field  Structure 

As  mentioned,  the  internal  fields  can  be  found  using  either  mode 
matching  or  a modal  expansion.  Since  the  numerical  results  are  found 
using  a computer,  it  is  advantageous  to  investigate  both  methods  for 
simplicity  and/or  speed  of  calculation. 

Within  the  cavity  the  fields  must  satisfy  the  boundary  conditions 
and  the  vector  wave  equation.  However,  one  major  difference  between 
modal  expansion  and  mode  matching  is  in  the  way  the  individual  mode 
terms  satisfy  these  conditions.  Although  both  methods  use  an  infinite 
set  of  orthogonal  functions  (modes),  a resonant  mode  expansion  requires 
that  each  mode  term  satisfies  the  boundary  conditions,  and  that  the 


total  field  satisfies  the  inhomogeneous  wave  equation.  Mode  matching 
requires  that  the  sum  satisfies  the  source  condition,  and  that  each  term 
obeys  the  source  free  wave  equation. 

The  vector  wave  equation  is  well  Known  to  be 


2-*  -►  V(V  • 7) 

V E + K E = jwuJ  - J 

1 


(2.2.1) 


where  is  the  complex  wave  number  in  the  dielectric.  For  the  thin 

cavity,  the  source  J has  some  special  properties.  For  the  microstrip 

♦ 

feed,  the  source  current  is  assumed  to  be  just  n X H,  where  H is  the  TEM 

H field  of  the  transmission  line.  For  the  coaxial  feed,  the  source  is 

taken  to  be  equal  to  the  current  density  on  the  center  conductor,  as  it 
existed  in  the  coax.  However,  this  source  is  simplified  by  considering 
the  current  distribution  to  lie  on  a cylinder  of  constant  radius.  Thus 
both  sources  are  considered  as  rectangular  distributions  in  $ and  as 

having  no  divergence.  For  simplicity  the  total  feed  current  is  assumed 

to  be  one  ampere,  flowing  into  the  port.  Both  source  distributions  are 
of  the  form  shown  in  Figure  2.3.  Therefore,  Equation  (2.2.1)  reduces  to 


2 

7 E + K E = jwpJ  . 
Z 1 z z 


(2.2.2) 


It  should  be  mentioned  that  for  a coaxial  feed,  no  attempt  has  been 

made  to  enforce  the  boundary  condition  that  E must  be  zero  on  the 

z 

coaxial  center  conductor.  Yet,  from  the  good  agreement  obtained  between 
the  calculated  and  measured  results  of  Chapter  3.  this  appears  to  be  of 


no  major  consequence.  This  could  be  attributed  to  the  fact  that  the 
conductor  length,  within  the  antenna,  is  very  short. 

A more  rigorous  treatment  might  be  to  consider  the  system  as  a 
cylindrical  cavity  with  an  offset  coaxial  feed.  The  internal  fields 
could  then  be  found  more  exactly  by  matching  modes  across  the  annulus  of 
the  feed.  However,  the  point  of  this  analysis  is  to  attempt  to  obtain 
approximate  analytic  results  with  a simple  but  perhaps  more  complete 
model  than  those  found  in  the  literature. 

As  in  most  problems  of  this  type,  the  Green's  functions  are  found 
first.  Hence,  the  feed  current  distribution  is  taken  as  a point  source 
and  the  resulting  fields  are  convolved  with  the  actual  source 
distribution  used. 

The  functional  form  of  the  source  free  resonant  modes  used  in  the 
modal  expansion  is  found  from 


2 2 

V ¥ + K ¥ = 0 . 

inn  mn  mn 


In  cylindrical  coordinates  this  gives 

I 


(2.2.3) 


¥ = J (K  p)cos(ni}>)  . (2.2.4) 

mn  n mn 


In  terms  of  these  modes,  the  electric  field  can  be  expanded  as  follows: 


H = — 
mn  top 


3f  „ at 
mn  mn  , . 

"T  ■ P - —7 — <f>  A 
p3<t>  3p  mn 


•] 


(2.2.6) 


Since  the  exterior  boundary  condition  is  that  H , = 0, 

ymn 


J ’ (K  a)  = 0 , 
n mn 


(2.2.7) 


where  a is  the  disk  radius.  Therefore,  the  resonant  frequencies  are 

those  where  (K  a)  is  a zero  of  J . 

mn  n 

Requiring  the  total  field  to  satisfy  the  inhomogeneous  wave  equation 
gives 


ju)pI5(*)6(p  - p’J/p'  = l A (K?  - O* 

* inn  I mn  mn  / 


m ,n 

=n 


where 


(2.2.8) 


? 

v y - -K  y « 
mn  mn  mn  ' 


(2.2.9) 


from  Equation  (2.2.3).  Here,  I is  the  source  current  in  amperes  and, 

from  the  symmetry  of  the  disk,  the  source  is  taken  to  be  at  (p* ,$  = 0)  . 

The  fact  that  the  {y  } are  orthogonal  can  be  used  as  an  advantage. 

1 mn 

Indeed  if 


C K J '(K  a)  + C J (K  a)  - 0 , 
j mn  n mn  2 n mn  ' 


(2.2. 10) 


then  from 
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/ J (K  p)J  (K  )pdp 
0 n mn  q pq 
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fj  ' OC  a)  J + (j  (K  a)  I 
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0 , otherwise  . 


(2.2.11) 


1 « •.  (m,n)=*(p,q) 

(K  a) 
mn 


we  can  say 


* T f f 6(«’)6(p  p‘) 

1 l i ’ 

P 4>  P pq 


y A (K?  - K2  ) f f V V pdpd* 
L mn  1 mn  1 ' mn  pq’ 


Performing  the  integrals  gives 


(2.2.12) 


2jupJ  (K  p*) 
n mn 


7ra2(l  + 6 )J  2(K  a)(K2  - K2  )|l  - — -«-} 

mO  n mn  1 mn  1 ...  .2' 
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Finally,  from  Equation  (2.2.4) 


(2.2.13) 


E ■ J A J (K  p)cos(n<|>)  . 
z u mn  n mn 


(2.2. 14) 


For  the  mode  matching  case,  one  needs  to  consider  the  homogeneous 
equation  in  each  source-free  region.  This  is 


12 


(2.2. 15) 


2 2 

7 1 + K y =0 

m 1 m 


where 


y = J (K  p)cos(m<*>)  or  y = Y (K  p)cos (m<t>) 
m m l m ni  i T 


(2.2. 16) 


or  any  linear  combination  of  Y and  J . The  source 

r.i  m 


J = If (4>) 6 (p  - p ' ) /p ' 


(2.2.17) 


can  be  restated  in  terms  of  its  Fourier  components.  Hence, 


t I5(p  - P1)  V u cos(mit’) 

J “ itp » 1 m (1  + 6 ) 

mo 


(2.2.18) 


m=0 
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Therefore,  the  source  current  can  be  considered  to  be  a sum  of 
cylindrical  current  sheets  located  at  p = p'. 

For  the  Green's  functions,  f(4>)  is  the  Dirac  delta.  Figure  2.4 
shows  the  two  regions  of  the  disk,  separated  by  the  cylindrical  current 
sheets.  Matching  of  the  fields  across  these  current  sheets  constitutes 
the  source  condition. 

In  region  1,  p'<p<a,  and 


oo 

E - y { A J (K  p)  + B Y (K  p)  } cos (m<(i) 
z j m m j m m i 

m=0 


(2.2. 19) 


At  p » a , 
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Figure  2.4.  The  disk,  showing  a source  current  distribution 
the  location  of  the  equivalent  cylindrical 
current  sheets,  and  the  two  source  free  regions 
of  the  mode-matching  analysis. 
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(2.2.20) 


With  the  orthogonality  of  the  cos(mij>),  the  result  from  region  1 is 


(A  J ' (K  a)  + B Y '(K, a)}  = 0 
m m 1 m m 1 


(2.2.21) 


where  the  derivative  is  with  respect  to  the  entire  argument. 

In  region  2,  Bessel  functions  of  the  second  kind  are  not  permissible 
since  they  require  a singularity  at  p = 0 . Hence, 


E = J C J (K  p)cos(m4>)  (2.2.22) 

Z2  m m j 

m=0 

Using  Equations  (2.2.19)  and  (2.2.22)  and  requiring  Ez  to  be  continuous 
at  the  source  gives 


A J (K  p')  + B Y (K  p')  = C J (K  p') 
m m j m m ] mm] 


where  the  orthogonality  of  cos  (m<}>)  has  been  used. 


At  the  source 


p X (H  - H )4>  = J 
1 2 


From  Equation  (2.2.6)  and  orthogonality, 


(2.2.23) 


(2.2.24) 


Equations  (2.2.21),  (2.2.23)  and  (2.2.25)  constitute  a determined 

set  of  equations  for  the  A , B and  C . Using  Cramer's  rule  for 

m m m 


determinants  yields 


A = - 
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jo)uY  '(K  a) J (K  p') 
tn  - j m - j 

2(1  +6  )J  ’(K  a) 
mo  m i 


jtopj  (K  p') 
m j 


2(1+6  ) 
mo 


juu  { J (K  p')Y  '(K  a)  - Y (K  p’)J  ' (K  a)  } 
m ] m j m j m j 


(2.2.26) 


2(1  +6  )J  '(K  a) 
mo  m i 

Finally,  Equations  (2.2.19)  and  (2.2.22)  can  be  used  to  find  the 
electric  field  anywhere  in  the  disk,  due  to  a point  source. 

The  results  from  the  resonant  mode  expansion  and  from  mode  matching 
can  be  checked  for  consistency  by  comparing  the  E field  Green's 
functions  from  the  two  methods.  From  the  theory  of  complex  variables, 
two  systems  are  equal  if  their  residues  are  equal  and  neither  system  has 
poles  at  infinity. 

Consider  any  point  within  region  2 of  Figure  2.4.  Multiplying 

2 2 

Equations  (2.2.14)  and  (2.2.22)  by  (K  - K ) and  taking  the  limit  as  K, 

1 mn  L 

approaches  K gives 
mn 


1 


I 


2 j coy J (K  p')J  (K  p) 
n mn  n mn 


a27ij  2(K  a) 
n mn 


[‘  ’ottf] 

» mn  -* 


2 2 

(K  - K ) 

• Um  1 M — 


K >K  J ’ (K  a) 
1 mn  n mn 


(2.2.27) 


J (K  p)J  (K  p’)Y  '(K  a)  . 


2 n mn  n mn  n mn 


But  since 


lim 


2 2 

(K  - K ) 2K 
1 mn 


mn 


K -+K  J ' (K  a)  aJ  "(K  a) 
j mn  n mn  n mn 


the  relationship  simplifies  to 


(2.2.28) 


tt(K  a) 
mn 


Y ' (K  a)  J (K  a) 
n mn  n mn 

J ' (K  a) 
n mn 


n 

1 2 
(K  a) 
mn 


(2.2.29) 


The  Wronskian  relation  that 


2 

it  (K  a) 
mn 


= Y ' (K  a)  J (K  a) 
n mn  n mn 


- J ' (K  a) Y (K  a)  , 
n mn  n mn 


(2.2.30) 


gives  the  requirement  that 


J "(K  a)  = J (K  a)  1 
n mn  n mn 


-> 


n 


(K 


mn 


(2.2.31) 


for  the  two  methods  to  be  consistent.  But,  Equation  (2.2.31)  is  true 

from  Bessel's  differential  equation,  and  inspection  of  Equations 

(2.2.14)  and  (2.2.22)  shows  that  no  poles  lie  at  K = °°.  Similar 

mn 

results  can  be  shown  for  any  point  in  region  1 of  Figure  2.4.  Thu3, 
mode  matching  and  a resonant  mode  expansion  yield  identical  results  for 
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the  simple  cavity  model  where  the  divergence  of  the  source  is  zero. 

Equations  (2.2.14)  and  (2.2.22)  indicate  that  evaluation  of  the 
internal  fields  requires  a double  summation  for  a modal  expansion  and  a 
single  summation  for  mode  matching.  Hence,  it  appears  that  less 
computational  effort  is  required  when  using  the  results  from  mode 
matching. 

The  fields  due  to  a source  with  a specific  feed  width  are  found  by 
convolving  the  Green's  functions  with  the  source  distribution.  Thus, 

P (<J>) 

E = * J F (p)  cos(m<f>) 

z 2a  u m 

m=0 

(2.2.32) 

r _ , , / .x  sin(ma) 

= > F (p)cos(m4>)  — 7— r — 

*•  m (.ma; 

m=0 


Here  P ($)  is  the  unit  "pulse"  of  width  2a,  and  F ($)  is  the  separated 
a m 

^'-dependent  function  found  in  equations  (2.2.19)  and  (2.2.22). 
Therefore,  the  electric  and  magnetic  fields  are  simply  the  Green’s 
functions  multiplied  by  {sin  (ma)}/(ma).  In  regions  1 and  2 the 
electric  fields  are  then 


<*>  sin(mi{>  ) 

E = y (A  J (K  p)  + B Y (K  p)}  cos  (m<p ) —7 — ~ 1 

zl  „ m m 1 m m 1 m4> 

m=0  w 


(2.2.33) 


E _ = / C J (K  p)  cos(mij)) 
z2  „ m m 1 


m=0 


sin(rvf>  ) 
\± 

m<() 


With  the  E field  known  throughout  the  disk,  the  time-averaged 


electric  stored  energy  can  be  found  from 


W = ~ f E dv  . 
e 2 

v 


(2.3.1) 


The  forms  for  the  fields  in  regions  1 and  2 of  Figure  2.4  are  stated  in 
Equation  (2.2.33)-  Since  the  fields  are  separable  in  <p,  d and  z,  the 
volume  integrals  are  easy  to  evaluate. 

Let  the  electric  field  be  represented  in  the  following  manner: 


Jzi  2 


l P (p)cos(miji)  , 


z2  2 L m 
m=() 


— Y Q (p)cos  (m>J>)  . 

m 


(2.3.2) 


Then,  the  electric  3tored  energy  is  expressible  as 


r C ( to  w ) 2 1 oo  co 

o r r r 


00  °°  b * a * 

} ) / Q Q pdp  + / P P pdp 

r.  « ,,  m n , m n 

m=0  n=0  0 p 


• / cos  (mi)))  cos  ( n<t> > d 4> 
0 


Since  the  last  integral  is  zero  except  when  m = n, 


t f (wp)2t71"<>  Ip' 


m=0  0 


W = -2“ l f | Q 1 2PdP  + / |P  |2pdp 


(2.3-3) 


(2.3-4) 


Each  integral  in  this  equation  can  be  evaluated  in  closed  form  since  the 
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and  are  just 


sin(m4>  ) 

P = {A  J (K  p)  + B Y (K  p)}  

m mml  mml  md> 

w 

sin(m<{>  ) 

Q = C J (K  p)  w 


(2.3.5) 


mm'  (1  + „)md) 

1 mO  w 


where  A ,B  and  C 
mm  m 

that 


are  given  by  Equation  (2.2.26).  It  is  well  known 


b 2 
/ U 2(K  p)pdp  = ^ 
m . 2 

a 1 


+ C1  - (Tjf]  K!(k1»>] 


(2.3.6) 


where  U is  any  linear  combination  of  j and  Y . By  takina  the 
m mm 

summation  inside  the  integrals  of  Equation  (2.3.4),  the  electric  stored 


energy  can  also  be  evaluated  using  a numerical  quadrature.  Thus, 


£ E TT  t ( (i)p  ) 2 


w = - 
e 


o r 


f l |Q  |2(1  + & )pdp 

0 m=0  m 


(2.3.7) 


a “J 

+/  £ |P  |?(l+6  )pdp 

I n m mo  l 

P m=0  J 

The  power  dissipated  in  the  dielectric  can  be  found  from 


P = o / | E | 2dv  . 
d 

But  this  is  Just 


P,  = 2o)6W  , 

d e ’ 

where  6 is  the  dielectric  loss  tangent. 


(2.3.8) 


(2.3.9) 
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The  copper  loss  is  approximated  by  initially  assuming  that  it  does 
not  change  the  internal  field  structure.  Therefore  it  is  calculated 


! 


using  the  internal  fields  derived  with  r "fectly  conducting  boundaries. 
Again,  this  method  is  plausible  since  the  copper  loss  is  very  small  with 
respect  to  the  reactive  power  within  the  disk.  The  copper  loss  is 
related  to  the  H field  tangent  to,  and  the  E field  normal  to  the  copper 
boundary . 

In  a good  conductor  the  skin  depth  is 


(2.3.10) 


where  o is  the  conductivity  in  mho-cm.  For  the  loss  due  to  the 
tangential  H field,  the  current  density  (J  ) within  the  copper  is 
assumed  such  that 


/ J dz  = J , 
„ v s’ 


(2.3.11) 


where  is  the  surface  current  density  that  would  exist  if  the 
conductivity  were  infinite.  is  the  unknown  volume  current  density 
flowing  parallel  to  the  surface,  with  form 


(2.3. 12) 


(2.3.13) 


so  that 

l r I ~Z/A 

'Js,e 

J = ~ h • (2.3.14) 

1 v A 

Since  the  fields  are  assumed  to  be  TEM, 

| J | = | H | . (2.3.15) 


The  power  dissipated  in  the  copper  is 


P = — / I J I 'dv 
cu,  a y v 


(2.3.16) 


so  that 


cu'  oA2  z=0 


o °°  , /.  a 2w 

— / e 2z/A  f / I J | 2 


r=0  *=0 


(2.3.17) 


But  this  is  just 


P = ■ - - W ? 

cu]  oAtii  h 


(2.3.18) 


where  the  factor  2 accounts  for  both  top  and  bottom  boundaries. 


* 
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The  copper  loss  due  to  the  E field  is  found  using  the  fact  that  the 


normal  D field  is  continuous  at 
permittivities  of  the  dielectric  and 

e.  = t e 

d or 

(2.3.19) 

e = e f 1 + -p—  1 , 

cu  o ^ Ju)LoJ 

and  for  a good  conductor 


the  copper  boundary.  Since  the 
copper  are 


jue 

o 


1 


) 


(2.3-20) 


then 


jwe  c 

E = — ^ E 
cu  o d 


(2.3.21) 


The  general  copper  loss  is 


Pcu  = 0 ; lEl  dv  • 

V 


(2.3.22) 


Since  the  E field  in  the  copper  varies  in  the  z direction  as 


E = e 

1 z cu 1 


-z/A 


(2.3.23) 


the  copper  loss  due  to  is 


CUo  o 


A(e  c io)?  ?v  a 

° r — / / | E , | ?dt 


(2.3.24) 


<p=0  p=0 


A 


where,  again,  the  factor  2 accounts  for  both  top  and  bottom  boundaries. 
This  is  expressible  as 


P 

cu2 


A2K  ;’W 

l_e 

otpA 


(2.3.25) 


The  total  power  loss  in  the  copper  is  now  seen  to  be 

2W,  + A?K  ?W 
hie 

p • 

cu  otpA 


(2.3.26) 


However,  the  skin  depth  is  so  small  that  the  loss  due  to  the  tangential 
H field  is  the  only  one  of  significance.  The  total  copper  loss  is  then 
approximated  by  Equation  (2.3*18). 

The  radiated  power  is  found  using  the  E field  at  the  edge  of  the 
disk.  From  Equation  (2.2.33),  the  field  at  the  edge  reduces  to 

. J (K  p ' )cos(m<t>)sin(m<t>  ) 

-Jti>u  ~ 

E (a , t> ) = 

z irK  a 

1 

where  $ is  the  azimuth  angle  with  respect  to  the  source.  Using  a 
Huygen's  surface  over  the  entire  cavity,  and  assuming  that  negligible 
currents  flow  on  the  exterior  surfaces  of  the  copper  boundaries, 


J ' (K  a)  (1  + 6 )m<J> 

m=0  mj  mo  w 


(2.3.27) 


I 


where  K 

s 


is  an  equivalent  magnetic  current  ribbon. 


(2.3.28) 


The  far  field  electric  vector  potential  is  well  known  to  be 


± _ - JV  t 2i  j(K  r-r') 

F ' f f Vie  ° ds* 

4,r  00  s 5 


(2.3.29) 


where  the  coordinate  system  is  as  in  Figure  2.2  and  the  factor  2 
accounts  for  the  ground  plane  image.  Assuming  the  thickness  to  be  small 
gives 


F = 


^ o 2it  jK  lx'sin(0)cos(4>)  + v sin(0)sin(<J>) 
e r t,  o 

— — / K e 
irr  Us 


ds',  (2.3.30) 


where  the  primed  variables  are  source  coordinates. 

The  vector  potential  can  be  found  as  the  superposition  of  its 
rectangular  components.  Since 


K = K a , 

S S ' 


(2.3.31) 


-xsin(j>’)  + ycos(<f>')  } 


(2.3.32) 


-jKor 

-Ituipe  ° °°  J (Kip ' ) sin(m<f>  ) 

^ m w 

2ir2rKi  m=0  J '(K,a)(l  + 6 )m<J> 

1 m 1 mo  w 


2n  „ jK  asin(0)cos(<fi,-(()) 

• f cos(m<J>)  {ycos(<J>' ) - xsin(<J>')le  ° d4>' 

0 


(2.3.33) 


«y 


Letting 


C(r)  = 


-JK  r 

-j2(DMte  ° 


Ati?K|  r 


and  a = KQa  sin(f>)  ^ 


(2. 3. 3*0 


and  using  the  fact  that 


ej„  ».(♦> , j (j)»J  (a)ej„ 

m y 

m=-°° 


(2.3.35) 


the  far  field  vector  potential  reduces  to 


; 


In  the  far  field 


E = V X F = - jK^r  X t + O^-^  f 


(2.3.39) 


so  that 


E = -jK^  - z sin(O)  ^cos 


s ( ()> ) F - s i n ( 4> ) F "1 

y *J 


(2.3.40) 


+ y cos(0)F  - x cos(0)F 
x v 


In  polar  coordinates,  the  E field  components  are 


E = cos(<f))F  - sin(<(>)F 
o y x 


E = 0 
z 


Ea  = -cos 


s(0)  £ cos(<f>)Fx  + sin(4>)Kv^j 


(2.3.41) 


Once  F^  and  F are  known,  the  far  electric  field  can  be  found  for  any  41 
and  0 . From  this,  the  radiation  pattern  can  be  calculated.  The 
radiated  power  is  found  from 


. n/2  2" 

P - — / / 

r no  0=0  $>=0 


/ |E|2r?  sin(e)ded<J> 


(2.3.42) 


where 


E 2 - E E + F.  E + E E_ 

xx  y y z z 


(2.3.43) 


jr 
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- . <*  .. 


The  integral  over  phi  can  be  carried  out 


in  closed 


form.  With 


orthogonality,  the  expression  for  the  radiated  power  can  be  further 
reduced  and  the  integral  over  9 can  be  numerically  evaluated. 


3.4  input  Impedance  and  the  Distributed  System 

In  general,  the  input  impedance  of  a rigorously  analyzed  system  can 
be  found  f rom 


P + j 2ai(W 

c 


V 


2 


W ) 
m 


(2.4.1) 


or 


P + j 2oi(U  - W ) 
e m 


(2.4.2) 


or  even  by 


Z 


in 


(2.4.3) 


Here,  P is  the  power  dissipated  in  the  system,  V is  the  driving  point 

voltage,  I is  the  source  current,  Z is  the  driving  point  impedance 

In 

with  no  other  ports  on  the  disk,  and  W and  W are  the  time-averaged 

e m 

electric  and  magnetic  stored  energies.  However,  for  the  disk  antenna, 
the  copper  loss  and  radiated  power  are  calculated  using  perturbation. 
Therefore,  the  above  equations  can  not  be  expected  to  agree  since  the 
input  voltage  is  insensitive  to  losses  calculated  in  this  way.  In  fact 


the  radiated  power  may  be  much  larger  than  any  other  loss,  so  that  the 


disagreement  between  the  above  equations  can  be  severe. 

An  extremely  efficient  method  for  calculating  the  input  impedance 
can  be  seen  by  considering  that  measurements  show  the  microstrip 
antennas  to  have  a high  Q.  Analytically,  the  quality  factor  of  any 
system  is  defined  as 


Stored  Energy 


Total  Power  Loss 


(2.4.14) 


Since  the  measured  Q is  very  high,  the  origin  of  the  power  loss  is 
unimportant.  This  indicates  that  the  losses  can  be  redistributed  with 
little  effect  on  the  internal  field  structure.  Since  the  input  voltage 
is  directly  related  to  the  E field  at  the  feed,  Equation  (2.4.3) 
indicates  that  the  input  impedance  will  be  relatively  unaffected  by  such 
a power  loss  redistribution.  This  result  may  be  paralleled  by  the 
simple  high  Q RLC  tank  circuit  whose  input  impedance  is  nearly 
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independent  of  where  the  loss  occurs. 

Once  the  copper,  dielectric,  and  radiated  power  losses  have  been 
calculated  as  in  Section  2.3,  they  can  be  distributed  throughout  the 
dielectric.  The  system  will  then  consist  of  a perfect  cavity  with  a 
modified  lossy  dielectric. 


The  losses  are  redistributed  by  modifying  the  dielectric  loss 
tangent.  From  Equation  (2.3.9),  the  modified  loss  tangent  is  defined  as 


The  complex  wave  number  for  the  modified  dielectric  is  then 


K1  eff  = “ \l(1  ~ j6eff)gE  (2-'*-6) 

and  the  E field  can  be  reexpanded  using  Equation  (2.2.33)*  The  input 
impedance  is  then  efficiently  evaluated  using  Equation  (2.4.3). 

The  input  voltage  can  be  found  using  the  E field  from  either  of 
regions  1 and  2,  and  the  fact  that 

V = E * zt  . (2.4.7) 


However,  since  the  voltage  is  averaged  over  the  feed  width  (due  to  the 
presence  of  the  feed  conductor),  the  final  expression  for  the  input 
voltage  is 


V 


in 


oo 
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m=0 


C J (K 
m m 
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P') 


sin*1(m4i  ) 
w 

(m$  J7 

w 


(2.4.8) 


It  is  interesting  to  note  that  the  quality  factor  of  the  disk  antenna  is 
just  the  reciprocal  of  the  effective  loss  tangent.  Indeed, 


Q 


ui(W  + W ) 
e m 


(2.4.9) 


However,  at  resonance.  Thus,  Equation  (2.4.5)  gives 


(2.4.10) 
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This  relationship  can  lead  to  some  major  computational  savings  since 
the  quality  factor  is  relatively  independent  of  the  feed  location.  Near 
the  m-th  resonance,  where  one  term  dominates  all  others,  it  is  clear 
that 


E(p=a)  <» 


J (Kip') 
m 

J '(K]a) 

m 


(2.4.11) 


But  the  radiated  power  is  proportional  to  the  square  of  the  electric 
field  at  the  disk  edge  so  that 


P a 
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J (KlP') 
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J ’(K,a) 
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(2.4.12) 


Also 
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e 


J (KlP') 
m 

J ' (K  j a) 
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(2.4.13) 


with  both  copper  and  dielectric  losses  being  proportional  to  Wg. 
Therefore,  as  long  as  the  resonant  term  dominates  all  others  in 
magnitude,  the  effective  loss  tangent  will  remain  nearly  constant.  This 
indicates  that  one  effective  loss  tangent  can  be  used  to  find  the  input 
impedance  for  any  feed  location  where  the  dominant  mode  is  excited.  The 
same  loss  tangent  can  be  used  for  any  frequency  near  resonance  since  one 
mode  is  still  dominant.  This  means  that  the  radiated  power  and  electric 
stored  energy  have  to  be  calculated  only  once  for  all  of  those  data 
points.  The  magnetic  stored  energy  is  never  needed  since  the  input 
impedance  is  equivalent  to  the  driving  point  voltage  when  the  source 
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current  Is  assumed  to  be  one  ampere. 


2,5  The  Disk  Antenna  with  Two  Ports 

The  general  multiple-port  disk  of  Figure  2.5  can  be  used  as  one  of 
many  elements  in  an  array  of  microstrip  antennas.  Another  use  for 
multiple  ports  is  in  exercising  some  form  of  mode  control  by  placing 
shorting  studs  (short  circuited  ports)  about  the  disk.  The  input 
impedance  of  a multiple-port  disk  can  be  approximated  by  using  the 
simple  model  of  the  previous  analysis.  However,  for  simplicity,  only 
the  two  port  disk  is  analyzed  here  since  the  analysis  readily 
generalizes  to  any  number  of  ports. 

A two-port  network  is  described  by 


1 

< 

1 

Zll 

Z l 7 

II 

CM 

> 

II 

z2, 

z22 

I2 

where  the  currents  are  directed  into  each  port.  The  open  circuit  Z 
parameters  have  been  chosen  since  all  are  calculated  with  one  of  the 
port  currents  equal  to  zero.  For  the  simple  cavity,  a port  with  no 
current  is  the  same  as  no  port  at  all,  so  that  each  parameter  reduces  to 
a one-port  calculation.  The  Z parameters  are 


2ij  - r | (2-5-2> 

1 I *0,  kilj 

**  A 

so  that  Z^n  is  equivalent  to  the  voltage  that  exists  at  the  location  of 
the  i-th  port  with  that  port  removed  (when  I is  taken  as  a one  ampere 


source).  From  equation  (2.2.33)*  the  voltage  is  known  for  any  location 
on  the  disk,  and  by  reciprocity,  = Z21  • f'r'om  Equation  (2.5.1),  the 
input  impedance  at  port  1 with  a loaded  second-port  is 


z .h..  , 

in  l1  "11  (Zl2  - Z22)  * 


where  Z^2  is  the  load  impedance  presented  to  the  second  port. 


(2.5.3) 


Lira? 


The  fields  in  the  disk  are  represented  by  an  infinite  series  as 
shown  in  the  analysis  of  Chapter  2.  It  has  been  found  that  the  use  of 
asymptotic  forms,  instead  of  truncation,  leads  to  more  accurate  analytic 
results  for  the  input  impedance.  Simple  asymptotic  forms  have  been 
found  for  both  the  input  voltage  and  the  electric  stored  energy. 
Radiation  due  to  the  non-resonant  modes  is  negligible  so  that  an 
asymptotic  form  need  not  be  found  for  the  radiated  power. 


The  advantage  of  using  the  asymptotic  forms  for  the  tail  of  a series 
is  that  a closed  form  expression,  or  one  which  converges  more  rapidly 


than  the  original  form,  can  often  be  found.  Using  Khummer's  well  known 


transformation,  any  series  S,  as  in 


Sl  = lim  y a 

p *“>  '•  m 

m=l 


(2.6.1) 


can  be  rewritten  as 


(2.6.2) 


where  the  and  are  the  actual  and  asymptotic  terms  respectively. 
As  indicated,  when  the  upper  limit  (P)  tends  to  infinity  the  sums  S1  and 
become  identical.  However,  by  using  Khummer's  transformation,  the 
rate  of  convergence  has  been  greatly  increased.  Therefore,  truncation 
of  S?  will  cause  less  error  than  the  truncation  of  S^. 

Investigation  of  the  equations  for  the  internal  fields  shows  that 
the  asymptotic  terms  must  use  the  large  order  approximations  for  the 
Bessel  functions.  For  large  m,  the  Bessel  functions  have  asymp- 
totic forms  of 


(2.6.3) 


From  Section  2.2,  the  electric  fields  in  regions  1 and  2 of  Figure 
2.4  are  known.  Substituting  the  asymptotic  forms  for  the  Bessel 
functions  gives 


In  simplifying  to  these  final  forms,  the  relations 
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were  used.  The  asymptotic  form  for  the  electric  stored  energy  is  found 
from  Equations  (2.6.5)  and  (2.2.33)  to  be 
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This  is  largest  when  p’=  a,  or 
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Even  for  this  largest  case,  4>,  is  always  large  enough  so  that  the  (W  ) 

w e m 

diminish  rapidly  with  increasing  m.  Therefore,  Equation  (2.6.1)  is  used 
to  calculate  the  electric  stored  energy.  Since  convergence  is  very 
rapid,  the  series  is  simply  truncated  at  some  order  N.  In  the  actual 
computer  program,  the  series  was  truncated  when  the  last  term  added  was 
less  than  .05?  of  the  accumulated  total. 

The  asymptotic  expressions  for  the  Z parameters  are  found  from 
Equation  (2.6.4).  As  pointed  out  in  Section  2.5,  the  Z parameters  all 
reduce  to  voltages  which,  in  turn,  reduce  to  electric  fields.  Averaging 
the  voltage  across  the  feed  width  gives  a Z ^ ( with  asymptotic  form 
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For  the  cross-port  parameters,  the  asymptotic  form  is  different  for  the 
two  regions  of  the  disk.  When  p i • p.  (where  the  j-th  port  is  the 
source  port), 
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Otherwise,  the  Pj  / P ^ term  is  inverted.  Inspection  of  these  equations 
shows  that  the  rate  of  convergence  decreases  as  p.  approaches  p. , and 
as  both  P^  and  Pj  approach  a. 

Rapidly  converging  forms  can  easily  be  found  for  these  asymptotic 

expressions  when  p = p = a or  when  p = ,,  . For  these  cases,  at 

i j i j 

least  one  term  in  the  brackets  of  Equation  (2.6.9)  is  identically  1. 
For  that  term,  the  sum  of  the  asymptotic  series  can  be  evaluated  quite 


rapidly,  using  Clausen's  integral.  Indeed, 
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But,  each  of  the  above  sums  is  known  to  be  the  integral  of  Clausen's 
integral,  found  in  any  major  handbook.  Numerical  evaluation  of  these 
sums  is  straightforward  and  is  outlined  in  Appendix  C. 

When  the  radial  feed  locations  are  not  equal,  the  series  of  (2.6.8) 
and  (2.6.9)  converge  at  least  as  fast  as  a geometric  series.  For  these 
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CHAPTER  3 

ANALYTIC  AND  EXPERIMENTAL  RESULTS 


3.1  Determination  of  6 and  o To  analyze  the  microstrip  disk,  the  copper 
conductivity  and  dielectric  loss  tangent  must  be  specified.  The 
conductivity  of  copper  is  well  known  to  be  580  Kmho-cm,  and  the 
dielectric  loss  tangent  is  given  by  the  manufacturer  of  the  dielectric 
in  the  form  of  a frequency  dependent  graph.  For  convenience,  it  is 
partially  reproduced  in  Figure  3.1.  However,  when  using  these  values  to 
calculate  the  input  impedance  locus,  th^. results  are  substantially  in 
error . 

To  illustrate  the  point,  a microstrip  disk  with  a radius  of  6.7  cm 
and  a thickness  of  1.5  mm  was  fed  at  a radius  of  3.36  cm,  and  analyzed. 
The  theory  indicates  the  lowest  resonance  to  be  near  810.61  MHz  since 
this  is  where  the  first  minimum  of  J ' (K^a)  occurs.  For  this  frequency 
the  loss  tangent  is  given  by  the  manufacturer  a3  0.00085.  The  Smith 
chart  locus  for  was  calculated  using  the  manufacturer's  loss  tangent 
and  the  conductivity  given  above.  For  comparison,  the  calculated  and 
measured  results  are  shown  in  Figure  3.2.  These  results  indicate  that 
more  loss  exists  in  the  system  than  the  theory  presently  accounts  for. 
At  this  point  it  was  decided  that  the  copper  loss  and  loss  tangent 
should  be  determined  experimentally. 

To  experimentally  determine  o and  6 , a bona-fide  rectangular 
cavity  was  constructed  by  closing  the  radiating  aperture  with  copper 
foil,  as  in  Figure  3-3.  With  the  bona-fide  cavity,  the  effective  loss 
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tangent  is 
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where  the  radiated  power  has  been  eliminated.  The  resonant  Q is  still 
given  by 


Since  P and  A are  given  at  resonance  by 
cu 
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Equation  (3.1.1)  leads  to 
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The  quality  factor  and  dielectric  thickness  can  be  measured,  leaving  6 
and  A (the  skin  depth)  as  two  unknowns.  By  constructing  two  similar 
cavities  with  different  thicknesses,  and  measuring  their  quality 
factors,  6 and  A can  be  determined  from  the  resulting  simultaneous 
equations.  Although  Equation  (3.1.3)  neglects  the  loss  in  the  copper 
foil  used  to  seal  the  cavity,  it  is  assumed  that  this  extra  loss  is 
negligible  since  the  cavity  thickness  is  much  smaller  than  any  other 


dimension . 


Once  the  skin  depth  is  determined  as  above,  the  conductivity  is 


found  from  the  relation 


(3-1.5) 


Unfortunately,  the  manufacturer’s  data  indicate  that  the  loss  tangent 
(for  the  Rexolite  2200  dielectric)  varies  with  frequency.  For  this 
reason,  6 and  o must  be  redetermined  when  the  frequency  is  changed 
substantially. 

The  rectangular  cavity  experiment  was  performed  for  a resonant 
frequency  near  428.00  MHz  and  the  results  are  shown  in  Figure  3.4.  The 
intersecting  lines  are  each  generated  by 


(3-1.6) 


With  the  materials  available,  cavity  thicknesses  of  1.5  mm  and  4.5  mm 
were  constructed.  Unfortunately,  these  thicknesses  give  a set  of 
equations  that  is  not  well  conditioned.  Because  of  this,  the  Q was 
measured  several  times.  The  locus  of  points  for  the  upper  and  lower 
limits  of  Q were  then  plotted.  In  Figure  3*1*,  the  resulting  region  of 
intersection  shows  the  uncertainty  of  <S  and  o . From  the  centroid  of  the 
shaded  area,  ,s  = .00083  and  A -9.2  um.  From  Equation  (3.1.5), 
o = 69.9  Kmho-cm.  This  method  can  be  performed  for  e.-ich  resonant  freq- 
uency . 

From  these  results  the  effective  conductivity  is  much  lower  than  the 
table  value.  This  discrepancy  may  be  due  to  surface  roughness  or 


possibly  to  a lossy  bonding  agent. 


3 • 2 Variation  of  Input  Impedance  with  Radial  Feed  Location 

From  the  simple  cavity  model  of  Chapter  2,  the  dominant  term 
contribution  to  the  resonant  input  resistance  can  be  predicted  to  vary 
as 
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where  (pq)  is  the  resonant  mode  index.  From  the  resonant  mode  expansion 
the  dominant  term  is 
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Since  the  dominant  term  is  the  only  one  that  contributes  a significant 
real  part,  equation  (3.2.1)  is  seen  to  be  true.  The  theory  was  tested 
for  disks  of  several  different  sizes,  and  the  agreement  between 
calculated  and  measured  results  was  very  good.  In  particular,  a disk 
with  a radius  of  6.7  cm  and  a dielectric  thickness  of  1.5  mm  was 
coaxially  fed  at  the  different  points  of  Figure  3*5.  For  the  analysis, 
a dielectric  loss  tangent  of  0.00135  and  a copper  conductivity  of  80.2 
Kmho-cm  were  used,  based  on  the  measured  cavity  results  of  Section  3.1. 
The  feed  width  was  taken  to  be  the  diameter  of  the  coaxial  center 
conductor.  Some  of  the  more  interesting  input  impedance  loci  are  shown 
in  Figure  3.6.  The  resonant  input  resistance  is  plotted  against  radial 


Figure  3.5.  Feed  location  used  to  determine  the  variatioi 
of  the  resonant  input  impedance  with  radial 
feed  location. 
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Figure  3.6.  Calculated  and  measured  first-mode  input 
impedance  loci  for  several  radial  feed 
locations  on  the  disk. 
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feed  location  for  all  eight  feed  points,  as  shown  in  Figure  3*7.  As 
seen,  the  calculated  and  measured  results  are  in  excellent  agreement. 

For  the  lowest  frequency  mode,  the  cavity  model  of  Chapter  2 has 
been  successful  at  predicting  both  the  resonant  input  resistance  and  the 
shift  of  the  entire  locus  above  the  real  axis  of  the  Smith  chart. 
However,  for  the  higher  frequency  resonances,  the  theory  predicts  too 
much  added  positive  reactance,  although  the  measured  and  calculated 
resonant  input  resistances  still  agree  very  well.  Some  results  for  the 
second  resonance  of  the  disk  of  Figure  3.5  are  shown  in  Figure  3.8.  The 
theory  predicts  too  Thuch  shift  above  the  real  axis,  regardless  of  the 
radial  feed  location.  In  Figure  3-9  the  resonant  input  resistance  is 
plotted  against  radial  feed  location  for  both  measured  and  calculated 
data.  This  Figure  shows  very  good  agreement  and  indicates  that  the 
losses  are  correctly  accounted  for.  For  a given  radial  feed  location, 
the  added  reactance  is  a function  of  the  feed  width,  as  Equation  (2.6.4) 
indicates.  By  increasing  the  feed  width  from  the  value  that  gave  good 
results  at  the  first  resonance,  the  results  from  the  higher  resonances 
can  be  made  to  agree  in  the  amount  of  locus  shift,  without  changing  the 
resonant  resistance. 

As  mentioned  in  Chapter  2,  no  attempt  has  been  made  to  satisfy  the 

boundary  condition  that  E be  zero  at  the  surface  of  the  coaxial  center 

z 

conductor.  It  may  be  that  this  has  contributed  more  error  to  the  higher 
frequencies  than  to  the  lower. 

The  results  for  the  higher  resonances  indicate  that  the  effective 
feed  width  of  the  coaxial  line  needs  investigation.  From  the  boundary 
condition  at  the  center  conductor,  the  electric  field  near  the  feed  has 
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Figure  3.9.  The  variation  of  second-mode  resonant  input 
resistance  with  radial  feed  location. 
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the  form  shown  in  Figure  3-10.  It  is  not  easy  to  see  how  this  field 
distortion  affects  the  input  impedance  contribution  of  each  mode.  For 
the  high-order  modes,  the  coaxially  fed  antenna  may  have  to  be 
considered  as  two  coupled  coaxial  cavities  as  suggested  in  Section  2.3. 


l-.l  Variation  of  input  imp^^nce_with  Shorting  stud  Location 

One  disadvantage  of  the  single-port  disk  antenna  is  that  the  feed 
must  be  taken  inside  the  disk  to  vary  the  input  impedance. 
Consequently,  the  feed  must  be  coaxial.  In  an  array  of  such  antennas, 
it  would  be  more  convenient  to  feed  each  antenna  with  a microstrip  line. 
This  would  eliminate  a lot  of  unwanted  coaxial  cabling  behind  the 
antennas . 

The  simple  cavity  model  predicts  that  a two-port  disk  antenna  can  be 
edge-fed  and  still  have  a variable  input  impedance  such  as  that 
discussed  in  Section  }.2.  From  Equation  (2.5.3),  if  the  second  port  is 
shorted  then 
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Considering  the  dominant  mode  contribution,  this  is 
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If  both  ports  are  at  the  same  radius,  then 


(3.3.3) 
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Thus,  the  resonant  input  resistance  is  seen  to  vary  with  radius  and  with 
separation  of  the  two  ports.  For  the  edge-fed  case, 

Z a J 2(K  a)  sin2(m<t>12)  , ( 3 . 3 • ^ ) 

in  p pq  > 

where  the  shorting  stud  is  also  on  the  edge. 

Equation  (3.3.4)  was  tested  with  a microstrip-fed  antenna,  with  a 
radius  of  12.68  cm.  The  first  mode  resonance  is  near  ^20  MHz.  Figure 
3.11  shows  the  disk  and  the  set  of  short  locations  used.  Some 
interesting  impedance  loci  are  shown  in  Figure  3.12.  Although  the 
calculated  shift  of  each  locus  is  incorrect,  the  resonant  input 
resistances  match  very  well.  Figure  3.13  shows  the  plot  of  maximum 
input  resistance  versus  separation  angle,  • The  agreement  is  shown 
to  be  very  good,  but  the  fact  that  the  impedance  loci  are  shifted 
incorrectly  indicates  that  the  feed  width  of  the  first  and  second  ports 
may  need  to  be  modified. 

Another  interesting  result  of  the  two-port  antenna  is  that 
additional  poles  are  introduced  into  the  system  where  Zf  2 = Z,,., , as 
Equation  (2.5.3)  indicates.  For  convenience  that  equation  is  repeated 
below. 
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When  Zj 2 i-3  a short  circuit,  a resonance  should  occur  near  where  Z 

passes  through  zero.  However,  for  an  edge  fed  disk,  Z ^ is  equal  to  Z^ 
when  the  feed  widths  are  the  same.  Note  that  the  "zero"  impedance  at 
port  2 is  converted  to  a "pole"  at  port  1.  This  was  verified  both 
experimentally  and  analytically.  For  the  disk  of  Section  3.2,  the  first 
resonance  of  Z22  occurs  at  810.61  MHz.  However,  as  shown  in  Figure 
3.14,  the  single-port  input  impedance  first  passes  near  zero  at  about 
860  MHz.  The  result  of  positioning  a short  circuited  port  at  $ = 45°  is 
also  shown  in  the  figure.  It  is  seen  from  the  figure  that  a second 
resonance  has  occurred  near  860  MHz,  as  predicted  by  the  theory.  That 
this  resonance  can  be  useful  is  shown  in  Section  3-5,  where  the 
radiation  pattern  is  investigated. 

3.4  A Circuit  Model  For  the  Disk  Input  Impedance 

In  the  analysis  of  the  disk,  it  was  shown  that  mode  matching  and  a 
resonant  mode  expansion  give  identical  results  for  the  simple  cavity 
model.  From  the  resonant  mode  expansion,  it  is  seen  that  the  input 
impedance  for  a one-port  disk  near  resonance  is 
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where  (pq)  is  the  resonant  mode  index.  Near  resonance,  the  first  term 
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in  the  above  equation  dominates  all  of  the  rest.  Considering  the 
dominant  term  alone,  the  input  admittance  has  a nearly  constant 
conductance  over  the  usable  bandwidth,  as  does  a parallel  RLC  tank 
circuit.  The  Smith  chart  locus  for  such  a contribution  is  shown  in 
Figure  3.15. 

The  analysis  of  Section  2.6  gave  an  asymptotic  form  for  the 
high-order  non-resonant  modes,  in  terms  of  the  mode  matching  results. 
From  Equation  (2.6. 1C)  and  any  handbook  on  Clausen's  integral,  it  is 
evident  that  the  sum  of  the  higher-order  non-resonant  asymptotic  terms 
contributes  a positive  reactance  to  the  input  impedance.  Experimental 
results  show  that  the  input  impedance  locus  is  nearly  a perfect  circle 
and  is  always  shifted  toward  inductive  reactance  for  the  resonances  that 
the  laboratory  equipment  is  able  to  examine.  The  result  would  seem  to 
be  that  the  sum  of  all  of  the  non-resonant  terms,  including  the 
lower-order  ones,  contributes  a positive  reactance  to  the  input 
impedance . 

Since  the  frequency  band  of  interest  is  very  narrow,  the  positive 
reactance  can  be  approximated  as  a linear  function  of  frequency,  and 
thus  as  an  inductive  reactance.  Adding  this  inductive  reactance  to  the 
dominant  mode  impedance  gives  the  shifted  locus  in  Figure  3.15. 

In  Sections  3*2  and  3 • 3 » it  was  shown  that  the  resonant  input 
resistance  of  the  dominant  mode  contributions  of  the  one-port  and 
two-port  disks  could  be  varied  by  changing  the  port  locations.  This 
behavior  is  quite  like  that  of  a tapped  RLC  tank  circuit,  where  the  tap 
location  is  changed.  Considering  these  results,  a simple  circuit  model 
for  the  one-port  and  two-port  disks  is  shown  in  Figure  3-16.  Equation 


Figure  3.15.  An  impedance  locus  for  an  RLC  tank  circuit, 
and  the  shifted  result  of  adding  the  effect 
of  non-resonant  terms. 
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(2.6.4)  indicates  that  the  series  inductance  in  Figure  3*16  is  a 
function  of  the  feed  width  and  radial  feed  location.  Since  the  sum  of 
the  non-resonant  contributions  has  a factor  w,  the  importance  of  these 
contributions  would  seem  to  be  a direct  function  of  frequency.  Indeed, 
the  loci  for  higher  resonances  are  shifted  above  the  real  axis  more  than 
those  for  the  lower  frequency  resonances  are.  However,  experimental 
results  indicate  that  the  shift  above  the  real  axis  is  less  than  the  u> 
factor  would  cause.  From  a term-by-term  inspection  of  the  input  . 

impedance,  some  terms  which  resonate  below  the  dominant  term  contribute 
negative  reactance. 

In  Figure  3.17,  the  points  of  maximum  resistance  are  shown  for  four 
of  the  radial  feed  locations  of  Section  3-2.  Although  the  points  seem 
to  lie  on  a line  of  constant  reactance,  Figure  3.18  shows  that  the 
reactance  increases  with  increasing  radial  feed  location,  as  Equation 
(2.6.4)  indicated.  The  results  are  somewhat  obscured  because  the  feed 
width  is  not  actually  constant  when  the  feed  is  moved  radially. 

From  this  simple  circuit  model,  "resonance"  of  the  antenna  can  be 
defined  as  the  point  of  maximum  input  resistance.  This  definition  allows 
the  antenna  and  the  RLC  tank  of  Figure  3-16  to  resonate  at  the  same 

■ 

frequency,  since  adding  reactance  to  the  RLC  tank  locus  does  not  change 
the  frequency  of  maximum  input  resistance,  as  seen  in  Figure  3-15. 


3.5  Radiation  Pattern  Investigation 

The  far  field  radiation  patterns  were  measured  and  compared  to 
calculated  results  for  the  one-port  disk  antenna.  The  patterns  for  the 
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multiple-port  disk  can  also  be  analyzed,  although  they  have  not  yet  been 
calculated.  The  radiation  patterns  are  calculated  from  equation 
(2.3.41),  where  the  E field  aperture  distribution  is  assumed  to  be  z 
directed  only. 

For  the  disk  of  Section  3.2,  the  first  resonance  is  near  810  MHz. 
The  calculated  and  measured  radiation  patterns  are  shown  in  Figures  3. 19 
to  3.24.  These  data  indicate  that  the  pattern  bandwidth  is  greater  than 
the  impedance  bandwidth.  Qualitatively,  the  calculated  and  measured 
results  match  quite  well.  However,  no  attempt  was  made  to  take  the 
effects  of  a finite  ground  plane  into  account.  For  frequencies  near  810 
MHz,  the  ground  plane  is  about  12  wavelengths  on  a side. 

For  a particular  mode,  the  radiation  pattern  is  substantially 
independent  of  the  size  and  resonant  frequency  of  the  disk,  as  long  as 
the  relative  permittivity  of  the  dielectric  is  held  constant.  With  this 
in  mind,  the  second  mode  pattern  was  investigated  using  the  disk  of 
Section  3-3-  For  this  disk,  the  second  resonance  is  near  710  MHz.  The 
calculated  and  measured  patterns  are  shown  in  Figures  3-25  and  3-26. 
Again,  the  agreement  is  seen  to  be  quite  good. 

The  measured  results  have  consistently  shown  higher  cross-polarized 
radiation  than  the  theory  accounts  for.  This  might  be  expected  since 
the  E field  at  the  edge  of  the  disk  is  not  exactly  z directed,  as  was 
assumed  with  the  simple  cavity  model. 

Although  no  radiation  patterns  were  calculated  for  a multiple-port 
disk,  some  were  measured  for  the  disk  of  Section  3-2.  A shorting  stud 
was  positioned  on  the  edge  of  the  disk,  45  degrees  from  the  microstrip 
feed.  The  angular  reference  for  all  two-port  measurements  is  taken  to 


Figure  3.19.  First-mode  radiation  patterns  for  the  6.7  cm 
disk,  f = 800  MHz,  <f>  = 0. 
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Figure  3.20.  First- mode  radiation  patterns  for  the  6.7  em 
disk,  f = 800  MHZ , <t>  = 00. 
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be  the  stud  location. 

The  radiation  pattern  was  measured  for  the  810  MHz  and  860  MHz 
resonances,  where  the  860  MHz  resonance  is  a result  of  the  zero-to-pole 
conversion,  as  described  in  Section  3.^.  The  810  MHz  pattern  results 
are  shown  in  Figures  3.27  and  3.28.  They  indicate  that  the  radiation 
pattern  is  "anchored"  to  the  shorting  stud.  This  suggests  that  the 

pattern  can  be  be  scanned  in  azimuth  by  repositioning  the  shorting  stud. 
However,  if  the  pattern  is  scanned  in  this  way,  the  input  impedance  will 
( vary,  as  describee1,  in  Section  3-2. 

The  radiation  patterns  for  the  860  MHz  resonance  are  shown  in 

Figures  3-29  to  3-31.  Figure  3-29  shows  the  azimuthal  radiation 

pattern,  and  indicates  that  the  pattern  maximum  is  aligned  with  the 

shorting  stud  (in  azimuth).  Figures  3.30  and  3-31  show  the  general 
shape  of  the  pattern  to  be  the  same  as  that  of  the  810  MHz  pattern. 
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CHAPTER  4 


SUMMARY 

In  this  thesis,  an  elegant  and  efficient  method  for  analyzing 
raicrostrip  antennas  with  simple  geometries  has  been  introduced.  An 
extensive  study  of  the  raicrostrip  disk  was  made,  and  analytic  and 
measured  results  were  compared  to  prove  the  validity  of  this  analysis 
method . 

A simple  but  accurate  analytic  model  for  the  microstrip  element  was 
developed  in  Chapter  2.  From  the  fact  that  the  radially  directed  disk 
surface  current  must  be  zero  at  the  edge,  and  from  the  close  proximity 
of  the  antenna  to  the  ground  plane,  the  antenna  element  was  modelled  as 
a closed  cavity  with  the  aperture  being  replaced  by  a perfectly 
conducting  magnetic  wall.  Since  the  antennas  are  very  thin,  the 
internal  fields  ware  taken  to  be  TEM  only.  With  the  assumption  that  a 
source  current  must  have  negligible  divergence  for  such  a thin  cavity, 
the  internal  fields  were  determined  and  confirmed  by  using  both  mode 
matching  and  a resonant  mode  expansion. 

An  extremely  efficient  method  for  calculating  the  input  impedance 
has  been  introduced  in  Chapter  2.  The  method  is  based  upon  the  fact 
that  the  input  impedance  of  a high  Q system  is  relatively  independent  of 
the  origin  of  any  power  loss.  Since  the  microstrip  antenna  has  a high 
Q,  the  radiated  power,  copper  loss,  and  dielectric  loss  can  be 
considered  as  a lumped  dielectric  loss,  without  substantially  changing 
the  calculated  input  impedance.  The  input  impedance  can  then  be  found 
by  finding  the  fields  inside  the  perfect  cavity  with  the  modified  lossy 


dielectric. 


The  microstrip  disk  exhibits  some  interesting  input  impedance 
properties,  as  the  analysis  above,  and  the  measured  data  of  Chapter  3 
can  verify.  Near  a resonance,  the  input  impedance  locus  is  an  almost 
perfect  circle,  slightly  shifted  above  the  real  axis  of  the  Smith  chart. 
Sections  2.3  and  2.6  show  that  the  diameter  of  the  impedance  locus  is  a 
function  of  the  dissipative  losses  and  the  radiated  power,  while  the 
amount  of  shift  above  the  real  axis  is  a function  of  the  effective  width 
of  the  source  current  distribution. 

With  the  simple  analysis  method  described  above,  it  can  be  shown 
that  the  disk  antenna  has  a very  useful  variable  input  impedance 
property,  specifically  useful  for  matching.  In  Section  3.2,  it  is  shown 
that  by  radially  varying  the  feed  location,  the  resonant  input  impedance 
can  be  varied  by  orders  of  magnitude.  In  Section  3-3,  it  is  shown  that 
a two-port  disk  can  have  the  same  property.  By  varying  the  angular 
separation  of  the  feed  point  and  a shorting  stud,  the  resonant  input 
impedance  can  again  be  varied  over  orders  of  magnitude. 

Section  3.1*  shows  that  the  input  impedance  for  the  one-port  disk 
antenna  can  be  described  very  accurately  and  elegantly  with  a simple 
circuit  model.  The  circuit  model  is  shown  to  consist  of  a series 
inductor  and  a tapped  RLC  tank  circuit,  where  the  RLC  tank  represents 
the  resonant  mode  contribution,  and  the  series  inductor  represents  the 
contribution  of  the  non-resonant  modes. 

Both  the  input  impedance  and  the  radiation  pattern  of  the  disk  were 
investigated  in  detail.  Chapter  3 shows  that  the  results  calculated 
from  the  simple  analysis  are  in  excellent  agreement  with  measured  data. 
However,  from  the  higher-order  mode  results,  it  appears  that  some 
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• j 

additional  work  must  be  done  toward  finding  the  effective  source  width 

j 

of  coaxial  feeds. 

The  accuracy  of  the  calculated  resonant  frequency  of  the  disk  was  , 

improved  somewhat  by  using  a method  repeated  from  the  literature  in 
Appendix  A.  With  the  method,  the  resonant  frequency  error  was  reduced 
from  more  than  2%  to  less  than  0.7%. 
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APPENDIX  A 


IMPROVEMENT  OF  CALCULATED  RESONANT  FREQUENCY 


From  the  cavity  model  of  Chapter  2,  the  resonant  frequencies  of  the 

raicrostrip  disk  are  found  to  be  those  frequencies  where  (K^a)  is  a zero 

of  J (K.a) . However,  since  this  is  a perfect  cavity  result,  the  field 
N i 

fringing  at  the  edge  of  the  disk  is  neglected.  In  reference  [4]  a 
method  has  been  described  to  approximate  the  effects  of  this  fringing. 
By  finding  the  dc  capacitance  of  the  disk,  an  effective  radius  can  be 
found,  and  is  always  larger  than  the  physical  radius.  This  larger 
radius  is  then  used  at  the  microwave  frequencies  for  all  calculations. 
From  the  above  reference, 


a = a • Ef 
ef  f 


(A  .1) 


where  a is  the  effective  radius,  and  EF  is  the  extension  factor,  and 
el  l 

is  defined  by 
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(A. 2) 


The  modified  resonant  frequency  is  easily  calculated  from 


f = 


V Z 
c o 


2na  ,,  e J 
ef  f r 


(A.  J) 





where  is  a particular  zero  of  Jn'(Z),  and  is  the  speed  of  light  in 
free  space.  As  an  example,  the  modified  resonant  frequency  is  compared 
to  the  original  for  the  disk  of  Figure  3-5.  The  measured  resonant 
frequency  is  found  from  Figure  3-6  as  792  MHz.  The  unmodified  resonant 
frequency  is  calculated  as  810.6  MHz.  This  is  in  error  by  2.3  %•  From 
equation  (A. 3),  the  modified  resonant  frequency  is  found  to  be  797.7 
MHz.  This  is  error  by  only  0.7  %.  Recall  that  from  the  definition  of 
resonance  in  Section  3.4,  and  from  Figure  3-6,  the  resonant  frequency  is 
independent  of  feed  location. 
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APPENDIX  B 

GENERATION  OF  COMPLEX  BESSEL  FUNCTIONS 

Numerical  techniques  which  are  used  to  generate  complex  Bessel 
functions  are  well-described  in  the  literature.  An  outline  is  presented 
here  for  completeness. 

Throughout  the  analysis  of  Chapter  2,  it  was  seen  that  the  wave 
number  (K()  was  complex.  Therefore,  the  Bessel  functions  having  in 
the  argument  were  also  complex.  The  Bessel  functions  are  found  using 
the  recursion  relations 


V ,(z)  + S'  X1 (z)  = - S'  (z) 
n-1  n+1  z n 


and 


S'  , (z)  - S'  . (z)  = 2 S'  ’ (z)  , 
n-1  n+1  n > 


(3.1  ) 


where  S'  is  either  Y or  J . 
m mm 

Since  the  large-order  asymptotic  expression  for  .1^  approaches  zero 

in  the  limit,  the  J are  found  by  using  a reverse  recursion  of  B.l.  The 
m 

recursion  is  initiated  at  some  high-order  N,  with  .J^  and  JN 
arbitrarily  chosen  as  (1+jO)  and  (O+jO),  respectively.  After  the 
recursion  is  completed,  the  requirement  that 

oo 

1 = J +2  l J (z)  (B.2) 

o , 2m 

m=l 

is  used  to  denormalize  all  of  the  J . Since  the  normalization  factor  is 

m 
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F = J +2 
o 


<■> 


J. 

m=l 


J9  (z) 

2m 


each  is  found  by  dividing  the  respective  recursion  result  by  the 
normalization  factor. 

Once  the  J are  found,  the  Y can  be  found  using  a forward  recursion 
of  B. 1 since  the  Y increase  with  increasing  order.  The  starting  values 
for  the  recursion  are 


Y (z) 
o 


k J2k(z) 

<->  -T- 


J 


and  from  the  Wronskian  relations, 


(B.4) 


ri(z)  = [Jl 


(z) 


Y (z)  - — 
O 11  z 


Vz) 


(B.5) 


Here,  the  principal  argument  is  used  for  ln(z)  and  care  is  taken  not  to 
specify  (z)  as  an  exact  zero  of  J^z),  or  of  JN1(z).  Also,  Y is  the 
constant  0.577 2. 

To  generate  the  J , an  order  must  be  chosen  at  which  to  start  the 
m 

recursion.  The  large-order  asymptotic  expression  for  is 


J (z) 
n 


(B.6) 


and  as  nw,  the  asymptote  approaches  zero.  The  initial  choice  of  (O+jO) 
for  JN  was  in  error  by  approximately 


88 


‘^In 


(B.7) 


since  J.,  was  assumed  to  be  zero.  By  specifying  an  acceptable  absolute 
N 


error  c , the  starting  order  can  be  found  by  numerically  solving  the 
equation 


N = i £ 


[1] 


N 


2 


(B.3) 


Four  digit  accuracy  was  desired,  so  • was  chosen  as  0.0001 
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APPENDIX  C 

EVALUATION  OF  CLAUSEN’S  INTEGRAL 


In  Chapter  2 it  became  neccesary  to  evaluate  the  sums 


f(0) 


i 

k=l 


sin  (k<>) 
k7 


(C.l) 


and 


g(U) 


«> 


); 

k=l 


cos  (ko) 

~"iTT 


(C.2) 


Equation  C.l  is  a restatement  of  Clausen's  integral,  found  in  any  major 
handbook,  and  is  just 


f(B)  = 0 


ln(0)  +2  l 
k=l 


C(2k) 

2k (2k  + 1) 


(C.  3) 


where  i,  is  the  Rieman-Zeta  function,  also  found  in  most  handbooks. 
Although  neither  C.l  nor  C.3  is  a closed-form  expression,  C.3  converges 
much  more  rapidly  then  does  C.l. 

Equation  C.2  is  related  to  C.l  since 


l 

k=  1 


cos  (kO) 
kT- 


"f  } d0  + c(3)  • 

0 k = l 


Thus  C.2  can  be  restated  as 


(C.  4 ) 
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